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ABSTRACT 

Development  of  a  three-component  long-period  borehole 
seismometer  of  small  diameter  (5  1/2  in.  O.D.)  is 
described.  The  feasibility  of  the  magnetic  suspension 
principle  has  been  demonstrated  and  several  instruments 
have  been  constructed  and  tested  in  the  laboratory. 
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I. 


INTRODUCTION 


The  present  report  covers  the  development  of  a  three-axis  long-period 
seismograph  suitable  for  operation  in  5-  to  6- inch  diameter  boreholes. 

This  program  relates  to  ARPA's  interest  in  further  improvements  of  the 
capability  to  discriminate  between  natural  and  man-made  tele^eismic 
events  of  small  magnitude  (m  ■  3.5  to  5.5  for  the  body  waves)  based  on 
the  comparison  of  their  apparent  short-period  and  long-period  magnitudes. 

While  the  detection  capability  provided  by  the  existing  short- 
period  borehole  seismometers  appears  to  be  quite  satisfactory,  the 
capability  of  the  long-period  instruments  cannot  be  fully  utilized  as 
long  as  they  are  Installed  in  shallow  vaults.  When  operated  near  the 
surface,  the  long-period,  horizontal-component  instruments  are  limited 
in  their  performance  by  the  high  level  of  ground  noise  at  low  frequencies. 
This  type  of  noise  is  believed  to  result  from  random  motions  of  ground 
blocks  generated  locally  by  environmental  factors  such  as  barometric 
pressure,  wind,  and  temperature.  In  order  to  avoid  this  situation,  the 
long-period  instruments  shoulo  be  placed  in  deep  boreholes  (500  feet  or 
deeper)  at  selected  sites.  Conventional  long-period  seismometers  are 
not  suited  for  installation  in  boreholes  because  of  their  large  size. 

Even  the  recently  developed  borehole  seismometers  require  boreholes  of 
large  diameters  (11  to  12  in.  I.D.)  which  tend  to  be  expensive. 

Development  of  the  magnetic  suspension  at  Arthur  D.  Little,  Inc. 
and  its  subsequent  applications  in  tiltmeters  and  other  instrument),  led 
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us  to  propose  the  development  of  borehole  seismometers  small  enough  to 
fit  in  boreholes  of  5-  to  6-inch  diameter. 

Practical  feasibility  of  such  instruments  has  been  demonstrated 
during  this  year’s  effort  and  the  instruments  have  successfully  under¬ 
gone  numerous  laboratory  tests. 

Even  though  it  was  not  anticipated  that  instruments  of  this  type 
would  necessarily  have  overall  performance  superior  to  that  of  high 
quality  station-type  instruments  which  are  in  existence,  it  appears 
that  the  performance  of  the  instruments  will  be  high  enough  tc  permit 
ground-noise  limited  operation  in  the  long  period  passband  at  most  of 
the  typical  sites. 
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II. _ SEISMOLOGICAL  BACKGROUND  FOR  THE  INSTRUMENT  SPECIFICATIONS 


The  present  program  relates  to  ARPA's  interest  in  further  improving 
the  capability  to  descr iminate  between  small  natural  and  man-made  tele- 
seismic  events  (  m  »  3.5  to  5.5)  by  comparison  of  their  magnitudes  at 
short  and  long  periods. 

While  arrays  using  existing  short  period  borehole  seismometers  are 
able  to  detect  everts  in  this  range,  long  period  horizontal  component 
instruments  are  usually  limited  in  their  performance  by  low  frequency 
tilting,  caused  by  local  environmental  factors  such  as  barometric 
pressure  changes,  wind  and  temperature  changes.  Much  of  this  tilt  noise, 
particularly  that  related  to  temperature  fluctuations  can  be  significantly 
reduced  by  careful  site  selection  and  by  placing  the  long  period  instru¬ 
ments  in  boreholes  several  hundred  feet  deep.  As  t'-e  cost  of  borehole 
preparation  climbs  rapidly  with  increasing  diameter,  practical  long-period 
instruments  should  combine  ruggedness,  reliability  and  small  (6  inches  or 
less)  diameter. 

The  operating  parameters  of  the  instrument  are  basically  determined 
by  the  comparative  power  distribution  in  the  spectra  of  the  signals  which 
one  is  trying  to  detect  and  of  the  microseismic  ground  noise.  The  sig¬ 
nals  being  sought  are  primarily  surface  waves  of  the  Love  and  Rayleigh 
types  generated  by  teleseismic  events.  Because  of  the  velocity 
characteristics  of  the  earth’s  crust  and  upper  mantle  surface  waves  from 
teleseismic  events  are  highly  dispersed,  wave  trains  of  many  minutes 
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duration  whose  energy  is  concentrated  predominantly  in  a  frequency  band 
between  approximately  0.025  and  0.10  Hz  (10  to  40  sec.).  Below  0.01  Hz, 
the  energy  of  the  signal  becomes  insignificant  except  for  the  largest 
natural  earthquakes  while  the  energy  of  the  ground  noise  begins  to  rise 
rapidly  from  its  minimum,  around  0.02  to  0.03  Hz. 

Larger  teleseisms  have  most  of  their  surface  wave  energy  concen¬ 
trated  near  the  microseismic  noise  minimum,  but  for  smaller  events  this 
energy  may  be  shifted  towards  somewhat  higher  frequencies.  We  show  in 
Figure  1  some  representative  data  collected  from  various  sources 
illustrating  typical  microseismic  power  spectra  in  the  range  0.01  to 
1.0  Hz.  While  it  is  customary  to  discuss  microseismic  noise  in  terms 
of  displacement,  for  consistency  with  the  discussion  in  following  parts 

of  this  report  we  have  expressed  the  data  in  terms  of  mean  square 

2  2 

acceleration  (cm/sec  )  per  Hz.  As  suggested  by  the  diagram,  the 
typical  microseismic  noise  spectrum  varies  greatly  from  place  to  place 
and  may  range  over  several  orders  of  magnitude  depending  on  meteoro¬ 
logical  conditions,  geological  conditions  and  distance  from  the  ocean. 

Although  most  of  the  noise  spectra  shown  in  Figure  1  are  based  on 
vertical  component  instruments,  on  the  high  frequency  side  of  the 
microseismic  minimum  the  magnitudes  of  the  horizontal  and  vertical  noise 
spectra  are  comparable  and  for  the  purposes  of  the  present  discussion, 
little  distinction  need  be  made  between  them.  At  long  periods,  however, 
the  similarity  is  not  so  pronounced.  Long  period  horizontal  seismometers 
are  notoriously  sensitive  to  tilts  whether  they  are  due  to  wind  loading, 
thermal  fluctuations,  barometric  pressure  changes,  or  earth  tides. 


II-2 


3rthur  2B.Htttlc.Jnr. 


FIGURE  1  COMPARISON  OF  ESTIMATED  ACCELERATION  POWER  SPECTRA  OF  SURFACE 
WAVES  FROM  NATURAL  TELESEISMIC  EVENTS  WITH  MICROSEISMIC  NOISE 
ESTIMATES  COMPILED  FROM  VARIOUS  SOURCES 


These  tilts  result  in  signals  which  by  themselves  are  indistinguishable 
from  long  period  surface  waves  and  their  frequently  large  amplitude 
imposes  stringent  requirements  on  the  linearity  and  dynamic  range  of 
the  instrument.  Although  the  relative  proportions  of  long  period  noise 
contributed  by  each  of  the  various  sources  is  still  a  subject  of  con¬ 
siderable  research,  it  is  clear,  both  from  our  own  work  and  that  of 
others,  that  extremely  small  temperature  fluctuations  on  the  pier,  or 
on  the  instrument  can  produce  major  effects.  By  placing  the  instruments 
in  deep  boreholes  these  temperature  effects  can  be  minimized  and  thus 
one  can  expect  to  remove  a  large  portion  of  the  locally  induced  noise 
leaving  only  the  signal  caused  by  propagating  seismic  waves. 

A  more  quantitative  comparison  of  the  relative  amplitudes  of  the 
teleseismic  signals  and  noise  is  shown  in  Figure  1  where  we  also  plot  the 
signal  level  from  earthquakes  of  various  magnitudes  M,  as  observed  on 
local  test  piers.  The  values  shown  were  obtained  by  first  estimating  the 
spectra  for  events  of  surface  wave  magnitude  around  7,  and  then  scaling 
to  smaller  magnitudes  assuming  that  the  shape  of  the  spectrum  has  not 
changed  in  the  region  of  interest.  From  this  diagram  it  is  clear  that 
except  for  the  most  ideal  sites  most  of  the  surface  wave  signals  from 
events  in  the  range  from  3.5  to  5.5  will  be  comparable  to  or  smaller 
than  the  microseismic  noise,  from  which  they  must  be  separated  by  various 
signal  processing  techniques.  We  also  show  in  Figure  1  some  preliminary 
data  on  the  microseismic  noise  at  two  local  sites  (Harvard  and  Bedford, 
Massachusetts)  which  we  have  been  considering  for  use  as  test  piers. 

At  both  these  sites  the  6-second  microseismic  noise  is  relatively  high, 
typical  of  locations  near  the  seacoast. 
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In  view  of  the  above  discussion  it  is  possible  to  reach  conclusions 
about  the  operating  characteristics  the  instrument  should  have.  0.01  Hz 
(100  sec.  period)  can  be  taken  as  a  lower  bound  for  the  operating  frequency 
band  of  the  seismometer.  At  the  upper  frequency  end  the  s.'gnal  tends  to 
overlap  the  large  microseismic  noise  peak  located  typically  between  0.15 
and  0.2  Hz.  Even  though  the  ground  noise  power  density  at  the  peak  may 
be  A  to  6  orders  of  magnitude  greater  than  at  the  minimum,  it  is 
desirable  (when  array  data  processing  is  considered)  that  the  instrument 
itself  have  the  capability  of  handling  this  large  dynamic  range  over  the 
entire  frequency  band  including  the  microseismic  peak.  Accordingly, 
the  upper  end  of  the  operating  frequency  band  may  be  placed  somewhere 
around  0.2  to  1  Hz. 

The  most  desirable  "center  frequency"  of  the  operating  frequency 
band,  usually  associated  with  the  natural  frequency  of  the  seismometer, 
is  determined  by  the  spectral  distribution  of  the  signal  and  the  ground 
noise;  but  in  practice  is  limited  by  the  requirements  of  instrument 
design  and  its  mode  of  operation.  We  may  estimate  that  the  center  fre¬ 
quency  will  be  most  likely  between  0.02  to  0.0A  Hz. 

The  absolute  sensitivity  of  the  seismometer  may  be  specified  in 
reference  to  the  minimum  signal  expected  and  the  level  of  the  micro¬ 
seismic  ground  noise  at  the  operating  location.  As  has  been  pointed 
out  above,  this  level  is,  of  course,  highly  variable  from  place  to  place 
and  depends  on  both  seasonal  and  occasional  meteorological  conditions. 
Depending  on  the  kind  of  intended  data  processing,  additional  margin 
between  ground  noise  and  instrument  noise  might  be  needed.  However,  as 
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a  first  estimate  it  appears  that  to  detect  surface  waves  from  events  of 

magnitude  Mg  =  3.5  to  5.5,  the  instrument  should  be  able  to  detect  sig- 

-12  2  2  / 

nals  of  average  spectral  density  of  10  (cm/sec  )  /Hz  over  a  bandwidth 
of  0.02  to  0.04  Hz.  This  is  approximately  the  level  of  the  lowest 
microseismic  noise  reported,  and  30  db  below  that  of  the  highest  values. 

Based  on  these  considerations  we  may  summarize  the  design  goals  of 
the  present  program  in  the  form  of  a  tentative  set  of  specifications 
as  follows: 


Type  of  Instrument: 


Operating  Frequency  Band: 
Ultimate  Sensitivity: 


Three-component,  6  in.  borehole 
seismometer  with  two  (identical) 
horizontal  component  modules  and 
one  vertical  component  module. 

0.01  to  1.0  Hz  (1  to  100  sec.  period). 

Approximately  10  mp  rms  ground  displace¬ 
ment  at  frequencies  between  .02  to  .04  Hz. 


Type  of  Suspension: 
Weight  of  Seismic  Mass: 

Damping : 

Readout: 

Mode  of  Operation: 
Dynamic  Range: 


Magnet ic . 

Approximately  1  and  10  gms  for  the  hori¬ 
zontal  and  vertical  modules,  respectively. 

Electromagnetic 

Displacement  type. 

Feedback. 

Approximately  80  db. 
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III. 


HORIZONTAL  SEISMOMETER 


A.  OPERATING  PRINCIPLES 


1_.  _ Diamagnetic  Suspension. 

The  horizontal  seismometer  employs  the  principle  of  diamagnetic 
levitation  to  suspend  a  seismic  mass  in  a  suitably  shaped  magnetic  field 
so  that  the  mass  is  constrained  to  move  in  the  direction  of  a  single 
horizontal  axis.  The  period  of  the  pendulous  motion  of  the  mass  is 
determined  essentially  by  the  depth  of  its  potential  minimum  due  to  the 
configuration  of  the  magnetic  field.  In  this  sense  the  elastic  force  of 
a  spring  or  the  gravity  force  acting  on  a  pendulum  is  replaced  by  a  mag¬ 
netic  field  force.  The  suspension  is  inherently  free  of  any  external 
(Coulomb)  or  internal  (anelastic)  friction.  However,  velocity- 
proportional  eddy  current  damping  can  be  introduced  at  will. 


The  field  configuration  utilized  in  the  horizontal  diamagnetic 
seismograph  is  characterized  by  a  strong  gradient  pointing  upwards, 
intermediate  gradient  in  the  lateral  horizontal  directions,  and  very 
weak  gradient  along  the  sensitive  axis  longitudinal  horizontal  directions. 
The  vertical  gradient  must  have  a  value  given  by  the  equation 


3H 

XvH  -r-2- 
y  3y 


mg 


(1) 


which  expresses  the  equality  of  the  magnetic  lifting  force  and  the  weight 
of  the  seismic  mass  (condition  of  levitation) ;  here  X  represents  the 
susceptibility  of  the  diamagnetic  seismic  mass,  v  the  volume  of  the  mass, 
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the  vertical  component  of  the  field,  and  mg  the  weight  of  the 

seismic  mass.  Since  X  has  a  small  numerical  value  (typically  about 

-6  x  10  ^  in  cgs  units  for  graphite,  which  is  the  best  available 

material),  the  value  of  H  (3H  /3y)  must  be  very  large,  typically  of  the 
8  2 

order  of  10  Oe  /cm.  Nevertheless,  gradients  of  this  magnitude  can  be 
readily  produced  with  commercially  available  permanent  magnets  in  work¬ 
ing  gaps  as  large  as  0.1  inch. 

The  lateral  gradients  are  obtained  by  shaping  the  field  in  a 
multipole  geometry,  or,  in  the  simplest  case,  by  an  edge-and-groove 
arrangement  of  magnet  polepieces  as  indicated  schematically  below: 


Fig.  2.  Polepiece  configuration  used  for  magnetic  suspension 
of  the  diamagnetic  seismic  mass  (cross  section). 
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The  sensitive  (horizontal)  axis  is  perpendicular  to  the  plane 
of  this  diagram  and  the  suspended  seismic  mass  (circular  in  cross 
section)  extends  some  forty  to  fifty  diameters  along  its  length.  The 
weak  longitudinal  gradient  along  the  horizontal  axis  is  obtained  by 
making  the  axial  (z)  component  of  the  field  initially  as  near  zero  as 
possible  and  then  introducing  a  small  perturbation  by  judicious  appli¬ 
cation  of  flux  shunts  or  disposition  of  magnets  along  the  polepieces. 

We  shall  outline  now  how  the  physical  parameters  of  the 
suspension  can  be  derived  from  the  properties  of  the  magnetic  field. 

The  field  configuration  according  to  Figure  2  is  not  readily  describable 
in  mathematical  form.  However,  fields  of  the  multipole  symmetry  which 
are  configurationally  similar  to  it  can  be  derived  from  a  magnetostatic 
potential  of  the  form 

n  «  C  J  (kr)  sin  n  0  cosh  (kz) .  (2) 

n 

This  potential  satisfies  Laplace's  equation  V2ft  =  0  and  therefore 
describes  a  physically  realizable  field.  Here  r  is  the  radial  coordi¬ 
nate,  6  the  azimuth,  and  n  an  integer  which  defines  the  periodicity  of 
the  field  about  the  z-axis  (n  -  number  of  pole-pairs  around  the  circum¬ 
ference);  Jn(kr)  is  the  nth-order  Bessel  function  of  the  first  kind,  and 
k  is  a  constant  which  determines  the  depth  of  the  potential  well  in  the 
z  direction.  The  field  components,  Hr>  H^  are  derived  from  equation  (2) 
as  -  3Q/3r  and  -r_1  3fl/30.  For  small  values  of  kr  and  kz,  equation  (2) 
has  the  approximate  form 

ft  *  A  r°  sin  n  0  (1  +  y  k^z^)  (3) 
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The  constant  A  can  be  expressed  in  terms  of  the  maximum  field  strength 


H  at  r  =  b  (the  distance  of  the  pole  tips  from  the  axis): 
m 


-  H  b(1"n) 
n  m 


The  potential  energy  V  of  the  diamagnetic  body  (a  cylinder  of  radius  a 

and  susceptibility  X)  suspended  in  the  field  is  determined  by  the 
2 

integral  of  H  over  the  volume  of  the  cylinder. 

The  potential  energy  is 


-  |  M  /  [nW-2  (i  +  i  „v) 


,  U  2n  2  2 

+  k  r  z  sin 


n  6  1 


dv 


(A) 


cyl . 


where  dv  is  an  element  of  volume  of  the  cylinder.  On  carrying  out  the 
integration  in  Eq.  (4)  for  the  case  where  the  cylinder  is  displaced,  an 
amount  h  radially  and  z  longitudinally,  one  can  calculate  the  longi¬ 
tudinal  spring  constant  K  from  the  formula 


K 


(5) 


One  can  show  that  for  n  =  2,  K  has  the  form: 

K  =  TT  £  (a/b) 2  I  x  |  H2  (h2  +  a2)  k2.  (6) 

Here  l  is  the  length  of  the  suspended  cylinder  and  H  is  the  maximum 

m 

field  at  the  polepiece  tips  (r  =  b) .  The  restoring  force  is  seen  to  be 

determined  by  the  magnitude  of  k,  but  is  also  dependent  on  H  .  There¬ 
in 

fore,  if  the  period  of  the  instrument  is  to  remain  constant,  the  product 
2 

|  x ,  H  should  not  vary  with  time,  temperature,  external  fields,  etc. 
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2. 


Readout  System  and  Instrument  Noise. 


The  readout  system  used  in  the  present  series  of  experimental 
seismometers  is  of  the  electrooptical  type.  The  quantity  sensed  by  the 
readout  is  the  displacement  of  the  seismic  mass  relative  to  the  frame  of 
the  instrument.  Because  of  the  proportionality  of  the  displacement  to 
the  inertial  force  acting  on  the  seismic  mass,  the  output  of  the  readout 
is  proportional  to  the  frame  (and  ground)  acceleration.  This  type  of 
readout  is  particularly  well  suited  for  the  diamagnetic  seismometer 
because  it  requires  no  bodily  contact  between  the  suspended  seismic  mass 
and  the  frame;  furthermore,  the  optical  readout  exerts  no  forces  on  the 
seismic  mass  in  the  direction  of  the  sensitive  axis. 

The  ultimate  limit  of  sensitivity  is  determined  by  the  "noise" 
inherent  to  the  instrument.  One  source  of  noise  is  the  photodetector 
used  in  the  optical  readout.  In  the  semiconducting  photodetectors  this 
noise  results  from  the  statistical  fluctuations  in  the  rate  of  generation 
and  recombination  of  the  current  carriers.  Unlike  Johnson  noise  its 
frequency  distribution  is  not  "white,"  and  its  magnitude  depends  on  the 
current  passing  through  the  device.  The  other  source  of  noise  is  that 
resulting  from  the  Brownian  motion  of  the  suspended  mass.  While  this 
noise  is  negligible  in  conventional  seismometers  employing  masses  of 
several  kilograms,  it  becomes  significant  in  the  present  seismometer  in 
which  the  seismic  mass  is  of  the  order  of  one  gram. 

It  can  be  shown  that  the  r.m.s.  value  of  the  "noise  displace¬ 
ment"  (z^)  of  a  mass  m  restrained  by  a  force  constant  K,  when  observed 
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over  a  frequency  range  from  0  to  “>  is  given  by  an  expression 


u 

where  =  (K/m)  is  the  undamped  natural  period  of  the  seismometer  and 
k^  =  1.38  x  10  ^  erg  sec  deg  ^  is  the  Boltzmann  constant.  In  order  to 
obtain  the  r.m.s.  value  of  the  noise  displacement  observable  in  a  given 
bandwidth,  it  is  necessary  to  consider  the  frequency  response  of  the 
seismometer. 


For  a  pendulous  system  of  a  mass  m,  natural  frequency  and 
damping  coefficient  a,  the  mean  square  noise  between  the  frequencies 
and  is 

fW2 

2  m  4kTa  i _ dw _ 

n  Timu^  J  1  +  2w^  ( 2a ^  -  1)  f  w^ 
o  w^ 


w  =  ii)/uo  is  the  dimensionless  frequency. 

This  result  can  be  derived  on  the  basis  of  an  equivalent  LCR 
circuit  in  which  the  inductance  represents  the  mass,  the  capacitance, 
the  restoring  force  and  the  resistance,  the  damping.  The  Johnson  noise 
in  the  damping  resistance  generates  a  charge  fluctuation  in  the  capaci¬ 
tance  which  corresponds  to  the  noise  displacement.  This  treatment 
yields  a  result  identical  with  that  obtained  from  the  kinetic  theory 
when  the  integral  is  taken  between  the  limits  0  to  °°.  Eq .  (8)  is 
therefore  a  generalization  of  Eq .  (7). 
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Evaluation  of  Eq .  (8)  for  a  seismic  mass  0.7  gms  under  typical 
operating  conditions  in  our  seismometer  gives  an  estimate  for  the  r.m.s. 

O 

value  of  the  Brownian  displacement  of  approximately  4.5  nm  (45  A).  We 

do  not  have  as  yet  typical  figures  for  the  noise  originating  in  the 

photodetector,  but  we  estimate  on  the  basis  of  preliminary  measurements 

that  it  is  smaller  than  the  Brownian  noise.  There  is  no  way  of  reducing 

that  part  of  the  photodetector  noise  which  lies  within  the  signal  pass- 

band.  The  Brownian  noise,  however,  can  be  reduced  by  reducing  the 

damping  of  the  viscous  or  resistive  type  acting  on  the  seismic  mass 

(z  -*  0  as  a  -*•  0  ,  see  Eq.  (8)).  Removal  of  damping  would  make  the 
n 

instrument  practically  useless  as  a  seismometer.  In  order  to  make  the 
instrument  perform  its  function  properly,  it  is  then  necessary  to  intro¬ 
duce  an  equivalent  of  damping  by  non-dissipative  means,  i.e.,  through 
the  use  of  a  feedback  circuit  which  controls  the  position  of  the  mass. 

3.  Feedback  Control. 

We  adopted  the  feedback  control  as  an  essential  part  of  the 
seismometer  originally  for  the  purpose  of  reducing  the  Brownian  noise 
as  explained  above.  On  closer  examination  one  finds  that  the  feedback 
control  offers  several  other  advantages.  First,  the  frequency  response 
characteristic  of  the  instrument  with  feedback  becomes  insensitive  to 
changes  in  the  period  of  the  diamagnetic  suspension.  Second,  the 
frequency  passband  of  the  feedback  instrument  can  be  varied  over  a  large 
range  by  varying  the  parameters  of  the  feedback  loop.  Indeed,  the 
response  can  be  made  flat  over  a  selected  frequency  range.  Finally,  the 
linearity  may  be  improved  and  the  dynamic  range  of  the  instrument  extended. 
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The  amount  of  electronic  noise  added  by  the  feedback  loop 
circuits  is  insignificant.  The  noise  originating  in  the  input  circuit 
(including  the  electrooptical  displacement  sensor  and  its  preamplifier) 
becomes  amplified  along  with  the  signal.  However,  it  can  be  shown  that 
the  signal-to-noise  ratio  remains  unchanged.  This  is  shown  in  a  detailed 
analysis  of  the  feedback  seismometer  in  Appendix  A.  The  principle  of 
operation  of  the  feedback-controlled  seismometer  may  be  described  as 
follows . 


Consider  a  seismic  mass  m  suspended  so  that  it  is  subject  to  a 
restoring  force  K'  (K  is  the  magnetic  spring  constant),  velocity 
proportional  damping  force  at,  and  an  external  (input  signal)  force  f(t), 
(Figure  3) . 


Fig.  3.  Schematic  diagram  of  the  horizontal  seismometer 
with  feedback. 
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4  *  z,  -  z  is  the  displacement  of  the  mass  relative  to  the  frame  of 
i  o 

the  instrument,  z^  is  the  coordinate  of  the  mass  in  absolute  (inertial) 
space,  z^  is  the  coordinate  of  the  frame  in  the  absolute  space. 

With  a  feedback  force  -f  =  AC  +  Be,  the  equation  of  motion  is 

me  +  (a  +  B)e  +  (K  +  A)e  =  f(t) .  (9) 

The  effect  of  feedback  is  to  change  the  spring  constant  and  the  damping. 

With  feedback  the  seismometer  will  have  a  pendulous  frequency 

^  (K+A) /m  and  damping  a+B.  If  A  and  B  are  independently  variable, 

the  frequency  response  and  damping  can  be  adjusted  at  will.  For  instance, 

the  so-called  maximally  flat  or  Butterworth  frequency  response  is 

2 

obtained  when  (a+A)  “2m  (K+B)  or,  if  A  and  B  are  large  compared  to 

2 

a  and  K  respectively,  A  =  2mB.  The  performance  of  the  feedback  instru¬ 
ment  is  then  essentially  independent  of  the  parameters  a  and  K  of  the 
original  seismometer. 

The  design  of  the  feedback  loop  is  most  readily  done  in  the 
frequency  domain.  The  basic  components  of  the  system  and  their  transfer 
functions  are  shown  in  Figure  4. 
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Optical 


Fig.  4.  Block  diagram  of  the  feedback  loop  circuit. 

Here  — 2 - - - 2  *s  t^e  acceleration-to-displacement  transfer 

s  +  2  6u>  s  +  w 
o  o 

function  of  the  seismometer,  is  the  displacement-to-voltage  constant 
of  the  optical  readout,  ^(s+a)  is  the  voltage  gain  of  the  amplifier 
over  the  frequency  range  of  interest,  and  is  the  voltage-to- 
acceleration  constant  of  the  feedback  element  or  forcer.  The  overall 

closed  loop  transfer  function  is 

v  K  K  (s+a) 

o  __ _  1  _  _  _ 

zi  s"  +  (2<S  ui  +  K  K  K.) s  +  u)2  +  K.K.K,  a  (10) 

o  1  2  J  o  12J 

When  K2  and  a  are  adjusted  to  appropriate  values,  the  frequency  response 

of  the  closed  loop  system  is  as  shown  in  Figure  5.  It  is  flat  to  a 

2 


frequency  u>  =  a  and  the  low  frequency  gain  is  —  + 

K3 


or 


KlK.a 
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FIGURE  5  EFFECT  OF  FEEDBACK  ON  THE  FREQUENCY  RESPONSE  OF  SEISMOMETER 


approximately 


K 


The 


is  also  shown  in  Figure 
this  case  is  — -r  . 


relative  frequency  response  of  the  seismometer 
5  for  comparison.  The  low  frequency  gain  in 


In  this  -espect  overall  feedback  serves  to  improve  the  linearity 
of  the  system  since  the  linearity  of  the  forcer  can  be  considerably  better 
than  that  of  the  seismometer  and  optical  readout. 


An  amplifier  transfer  function  of  the  form  K(s+a)  is  not 
realistic  since  the  gain  cannot  increase  with  frequency  without  limit. 
Therefore,  in  practice,  the  amplifier  transfer  function  will  be  of  the 
form  K(s+a) / (s+b) .  Analysis  of  a  feedback  with  this  transfer  function 
is  somewhat  more  complicated  than  the  previous  case,  but  the  important 
result  one  obtains  is  that  damping  of  70%  or  better  can  be  obtained  over 
a  range  of  gain  which  increases  with  the  frequency  ratio  b/a.  This 
result  is  shown  in  Figure  6. 


Fig.  b.  Gain  range  required  to  obtain  damping  greater  than  70% 
as  function  of  gain  step  ratio. 
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B.  INSTRUMENT  CONSTRUCTION 


The  critical  component  of  the  magnetic  suspension  for  the  long- 
period  horizontal  seismometer  is  the  set  of  polepieces  that  determine 
the  field  according  to  principles  outlined  in  Section  III,A1.  The 
easier  part  of  the  problem  is  the  cross-sectional  geometry  of  the  pole- 
pieces  as  outlined  in  Figure  2.  This  was  worked  out  earlier  and 
optimized  by  means  of  an  analog  field  plotter  so  that  a  maximum  support¬ 
ing  force  was  obtained  for  diamagnetic  cylinders  of  the  selected 
diameter  (0.086  in.)  with  practically  available  permanent  magnets. 
Shaping  of  the  field  along  the  sensitive  axis  is  far  more  difficult 
since  for  long-period  instruments  the  required  field  variation  from 
center  to  end  is  typically  one  part  in  ten  thousand  or  less.  The 
initial  approach  was  to  give  the  gap  a  slight  curvature  along  the  axis 
by  machining  the  edge  of  the  lower  polepiece.  As  it  proved  difficult 
to  increase  the  period  in  this  way  beyond  approximately  six  seconds, 
we  decided  to  make  the  polepieces  perfectly  straight  and  as  long  as 
compatible  with  the  overall  size  of  the  instrument.  By  making  use  of 
the  lengthwise  flux  distribution  in  the  polepieces  and  by  adjusting  the 
length  of  the  graphite  rod,  we  found  it  possible  to  make  stable  suspen¬ 
sions  with  a  natural  period  of  the  seismic  mass  from  16  to  30  seconds. 

The  homogeneity  of  the  polepiece  material  for  long-period 
seismographs  turned  out  to  be  rather  critical  for  reasons  of  imperfec¬ 
tions  present  in  common  stock  of  soft  magnetic  steels.  Cold  rolled 
Armco  iron  carefully  inspected  for  flaws  and  voids  was  found  adequate. 
The  polepieces  are  fabricated  from  a  round  stock  in  one  operation 
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without  re-chucking  while  held  between  a  dividing  head  and  a  center 
throughout  the  entire  machining  operation.  The  polepiece-and-magnet 
assembly  was  shown  in  Figure  2  of  our  Second  Quarterly  Report  dated 
31  December  1967. 

Overall  view  of  the  complete  horizontal  seismometer  is  shown  in 
Figure  7.  The  upper  part  is  the  seismometer  proper(S)  ,  the  middle  part 
is  the  remote  levelling  jack  screw  (RL) ,  and  the  bottom  part  is  the 
piezoelectric  forcer  or  tilt  platform  (PZT) .  The  instrument  is  housed 
inside  a  4.25  in.  I.D./4.75  in.  O.D.  iron  casing  (C)  with  gas-tight  lids 
and  cable  feedthrough.  Besides  its  protective  function  the  casing 
serves  as  a  magnetic  shield.  Assuming  permeability  u  =  1000  its  shield¬ 
ing  factor  should  be  better  than  50. 

The  "C"  shaped  Alnico  V  magnets  used  in  the  seismometer  are 
2.50  in.  wide  and  are  shaped  to  fit  the  angular  polepieces  of  3.25  in. 
axial  length.  The  "C"  magnets  adhere  strongly  to  the  polepiece 
structure  by  magnetic  attraction  and  the  assembly  is  fastened  to  the 
base  by  means  cf  clamps  and  set  screws  on  the  magnets. 

The  graphite  rods  are  machined  from  0.125  in.  diameter  spectro¬ 
scopic  grade  graphite  to  a  diameter  of  0.086  in.  with  a  diamond  tool. 

The  rods  must  be  thoroughly  cleaned  by  boiling  in  1:1  hydrochloric 
acid,  ultrasonic  cleaning  in  detergent,  and  boiling  in  several  changes 
of  distilled  water.  The  graphite  easily  picks  up  enough  iron  contamina¬ 
tion  from  the  machining  operation  to  make  it  behave  erratically  in  the 
magnetic  field. 
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FIGURE  7  VIEW  OF  THE  COMPLETE  HORIZONTAL  SEISMOMETER 


The  optical  readout  (Figure  8)  consists  of  a  light  source  (I.) ,  a 
light  obstructing  "flag"  (F)  attached  to  the  graphite  rod  and  dual 
photodetector  (PD)  of  the  cadmium  sulfide  (CdS)  photoresistive  type. 

The  lan.p  is  a  selected  T-l  size,  5  volt  bulb  rated  for  a  10^  hours' 
life  service.  The  CdS  photodetector  is  a  standard  CL  705  L/2  Clairex 
cell.  The  lamp  and  cell  are  provided  with  glass  light  pipes.  The 
light,  after  passing  past  the  flag,  is  conducted  to  the  two  halves  of 
the  photoconductive  cell  by  curved  (90°)  rectangular  light  pipes  (LP) . 

This  arrangement  permits  us  to  place  the  photocell  proper  inside  one 
of  the  magnets  Co  save  space. 

Levelling  of  the  instrument  after  it  is  placed  in  a  borehole  is 
provided  by  a  motor-operated  jack  screw  assembly  located  immediately 
below  the  seismometer  assembly.  The  motor  leveller  has  a  range  of 
approximately  ±2°  and  operates  on  a  one-  to  five-volt  D.C.  bipolar  source. 
The  motor  is  a  PM-field  D.C.  motor  with  a  560:1  gear  reduction.  There  is 
a  20:1  worm  gear  reduction  following  this  which  turns  a  40-pitch  pre¬ 
cision  ground  jack  screw.  With  five  volts  applied,  the  jack  screw  moves 
about  .001  inches  per  minute.  This  corresponds  to  a  tilt  of  400  micro¬ 
radians  per  minute.  We  find  that  the  motor  leveller  is  sensitive  enough 
in  its  operation  to  be  used  as  an  intermediate  levelling  device  such  as 
a  drift  corrector  to  keep  the  trace  on  the  chart  paper.  This  is  presently 
accomplished  by  using  high  and  low  limit  switches  on  the  recorder  which 
activate  the  motor. 
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FIGURE  8  VIEW  OF  THE  OPTICAL  READOUT  OF  THE  HORIZONTAL  SEISMOMETER 
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The  piezoelectric  tilt  table  located  under  the  motor-operated 
leveller  is  used  in  the  feedback  loop  to  deliver  the  controlled  force 
acting  against  the  inertial  force  from  the  input  acceleration.  If  the 
table  is  tilted  by  an  angle  0,  the  force  acting  on  the  seismic  mass  is 
mg  sin  0,  or  approximately  mg  0.  The  tilting  is  obtained  by  three 
stacks  of  piezoelectric  ceramic  discs  arranged  so  that  application  of  a 
+100  volt  control  voltage  tilts  the  seismometer  by  +5  microradians.  In 

terms  of  the  feedback  circuit.  Figure  4,  this  corresponds  to  a  forcer 

-5  -2  -1 

constant  ~  5  x  10  cm  sec  volt 

A  typical  horizontal  seismometer  has  a  natural  frequency, 

'  0.4  rad/sec  (16-second  period)  and  is  lightly  damped  (6  :  1/8),  so 
that  its  transfer  function  is 


z 

o 


_ 1 _ 

s^  +  .08s  +  .16 


(ID 


The  optical  readout  constant,  ,  is  about  20  volts/cm.  The 

feedback  amplifier  was  designed  with  both  the  gain  and  the  frequency  of 
the  gain  step  adjustable  over  wide  ranges.  At  present  the  gain  step  is 
fixed  at  10  (20db)  to  allow  for  a  2:1  gain  variation  with  acceptable 
damping.  The  operating  settings  result  in  a  voltage  transfer  function 
for  the  amplifier  approximately 


v 

o 


zi 


10 


4  s  +  2 

s  +  2  o' 


and  the  closed  loop  instrument  is  slightly  underdamped. 


The  circuit  diagram  of  the  horizontal  feedback  circuit  is  shown 
in  Figure  9. 
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OUTPUT  TO  PZT 
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FIGURE  9  HORIZONTAL  SEISMOMETER  FEEDBACK  CIRCUIT 


IV.  VERTICAL  SEISMOMETER 


A.  OPERATING  PRINCIPLE 

The  vertical  seismometer  makes  use  of  a  magnetic  field  suspension 
of  a  type  quite  different  from  that  employed  in  the  horizontal  units. 
Analysis  shows  that  the  vertical  suspension  of  the  diamagnetic  seismic 
mass  having  the  required  long  period  would  require  extremely  close 
tolerances  on  the  exact  shape  of  the  field  and  it  would  support  still 
smaller  masses  than  those  used  in  the  horizontal  diamagnetic  suspension. 

The  suspension  chosen  for  the  vertical  seismometer  is  of  a  hybrid 
type  consisting  of  magnetic  dipole  suspension  for  the  seismic  mass  and 
a  pair  of  diamagnetic  guides,  schematically  outlined  in  Figure  10. 


Fig.  10.  Schematic  diagram  of  the  vertical  magnetic  suspension. 
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The  seismic  mass  contains  a  permanent  magnet  of  a  dipole  strength  p 
which  is  suspended  in  a  magnetic  field  having  a  vertical  gradient  such 
that  the  upward  force  acting  on  the  dipole  is  equal  to  the  weight  of 
the  seismic  mass: 

3H 

P  =  mg.  (12) 

The  suspension  is  unstable  inasmuch  as  the  dipole  has  the  tendency  to 
flip  over  at  a  slightest  departure  from  exact  symmetry.  This  is  pre¬ 
vented  by  the  two  diamagnetic  guides  which  act  as  noncontacting  (fric¬ 
tionless)  sleeve  bearings.  They  consist  of  graphite  cylinders  supported 
by  magnetic  fields  of  a  quadrupole  configuration. 

The  entire  structure  is  symmetric  about  a  vertical  axis  and  the 
seismic  mass  is  constrained  to  move  along  this  axis  in  the  vertical 
direction. 

The  main  field  supporting  the  mass  may  be  described  by  a  magnetic 
potential  fi  represented  by  a  sum  of  cylindrical  narmonics  satisfying 
the  Laplace's  equation.  If  we  consider  only  the  first  three  harmonics, 
the  magnetic  potential  P  has  the  general  form: 

P  *  Az  +  B(2z2  -  r2)  +  C(2z^  -  3r2z)  (13) 


where  A,  B,  and  C  are  adjustable  constants. 

The  field  components  are: 

H  -  -  -7^  ■  -  A  -  4Bz  -  C(6z2  -  3r2) 

z  6  z 

H  *>  — m  2Br  +  6Crz 

r  6r 
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The  potential  energy  of  the  dipole  in  the  field  H  is 


-(y.fi)  =  cosQ  +  sin9  cos4^  (16) 

where  0  is  the  inclination  of  the  dipole  axis  to  the  z-axis  and  4  is  the 
azimuth  of  the  inclined  dipole.  In  (16)  the  north  pole  of  the  dipole  is 
assumed  to  be  directed  downwards. 

The  potential  energy  of  the  two  graphite  cylinders  is 

V'  =  x'v'H'2  (17) 

where  x'  is  the  diamagnetic  susceptibility  of  graphite,  v'  is  the  total 

volume  of  the  two  cylinders,  and  H’  is  the  quadrupole  field  in  each 

bearing.  The  magnetic  potential  for  these  quadrupole  fields  is 

2 

fl  ®  r  cos  2<(i  and  the  field  components  are  *  r  cos  24,  «  r  sin  24. 

Thus  H'  *  H'  r'/a  (18) 

m 

where  is  the  magnetic  field  at  the  quadrupole  pole-faces  of  radius  a. 

The  displacements  of  the  upper  and  lower  diamagnetic  cylinders  due 
to  the  tilt  of  the  suspended  assembly  are  given  by 

*=  r2  +  2£r  sin0  cos4  +  £2  sin20  (19) 

r|2  =  r2  -  2£r  sin0  cos4  +  £2  sin23  (20) 

where  r^,,  r^  are  the  radial  positions  of  the  centers  of  the  two  cylinders, 
r  is  the  position  of  the  permanent  dipole,  and  21  is  the  distance 


IV- 3 


3rtbur  Sl.^iittlr.JInr. 


between  the  centers  of  the  cylinders.  The  potential  energy  of  the  two 
cylinders  is  therefore,  from  (17)  and  (18\ 


V'  - 


x'v'h;2  o  o 

- —  (r’2  +  r'2) 

2  2  1  ' 


2x'v'H'2  2  2  2 

- - -  (r  +  l  sin  0) 


(21) 


Thus,  the  total  magnetic  potential  energy  of  the  suspended  system  is 


2x'v'H'2 

V  -  - r-21-  (r2  +  l2  sin20)  +  u  (H  cos0  +  H  sin0  cos^)  (22) 

2  z  r 


When  the  supporting  field  is  defined  by  the  potential  of  Eq.  (13),  the 

components  H  and  H  are  given  by  (14)  and  (15).  Therefore,  V  becomes 
z  r 


2XVH'2  9  9  9 

V  _  - (r  +  l  sin  0) 


+  ucos0  A  -  4Bz  -  C(6z2  -  3r2)J 


+  usin0  cos4>  [2Br  +  6Cr^ 


(23) 


We  can  determine  the  constants  A,  B,  and  C  by  considering  the  case 
r  *  0,  0*0.  Then 


V 


uA  -  4uBz  -  6uCz 


2 


(24) 
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On  applying  the  conditions  at  z  ■  0,  that 


H  - 

Ho 

(25) 

9  V 

-  ts 

9z 

w 

(26) 

2 

9  V 

-  '  S 

K 

(27) 

9z2 


we  obtain 


-  H0  ,  B 


W.  c 

4u  ’  L 


K _ 

12u 


(28) 


The  potential  energy  of  the  suspended  sy  tem  can  therefore  be 
written  as: 


2  X  ' v ' H '  ,  2  2 

v  .  - _JL.  (r^  +  ^  sin^e) 

a 

+  cosG  £  UHq  "  Wz  +  -  3r2)  J 

+  sin6  cos$  (j  r  -  y  rz)  (29) 

For  small  values  of  the  variables  z,  r,  6,  and  $  we  find 


(30) 
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Note  that  $  enters  only  in  terms  higher  than  second  order. 


Stability  requires  that  V  must  increase  from  its  equilibrium 
value  of  uHq  for  any  arbitrary  combination  of  z,  r,  and  0.  The  neces¬ 
sary  conditions  are 


2x'v'H'2 

m 

2 

a 


> 


K 

4 


2x'v'£2H'2  uHa 

_ m_  _ y 

2  2 
a 


(31) 


(32) 


(33) 


Condition  (31)  has  to  do  with  pure  radial  displacement,  (32)  with 
pure  rotation,  and  (33)  with  combined  rotation  and  radial  displacement. 


It  is  to  be  noted  from  (30)  that  although  there  is  coupling 
between  r  and  0,  there  is  no  interaction  of  these  coordinates  with  z. 

This  fortunate  state  of  affairs  is  only  realizable  with  the  field  defined 
by  (13).  Any  field  more  complicated  than  this  will  produce  cross¬ 
coupling.  On  substituting  from  (28)  into  (14)  we  find  for  the  required 
field  on  the  axis: 


H 

z 


0  u  2u 


(34) 
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A  field  of  this  type  has  equipotential  surfaces  of  the  form  of  a  hyper¬ 
boloid  of  revolution  with  asymptotes  defined  by 

z  =  t  — —  r.  (35) 

JT 

This  is  tne  field  we  actually  use  in  the  vertical  seismometer. 

A  field  of  this  geometry  can  be  realized  with  a  set  of  polepieces 
having  equipotential  contours  according  to  Eq.  (34).  The  polepiece  set 
having  a  shape  of  a  hyperboloid  of  revolution  is  energized  by  permanent 
magnet  rings  placed  at  appropriate  locations  around  the  working  gap. 

Even  though  this  type  of  suspension  can  suspend  considerably  larger 
seismic  masses  than  the  diamagnetic  suspension,  we  followed  the  same 
general  philosophy  in  the  overall  design  of  the  vertical  seismometer. 

In  particular,  we  decided  to  provide  the  instrument  with  feedback  con¬ 
trol  for  the  same  reasons  as  with  the  horizontal  instrument.  Detailed 
analysis  of  the  vertical  feedback  seismometer  is  given  in  Appendix  A. 

B.  INSTRUMENT  CONSTRUCTION 


The  vertical  magnetic  suspension  system  forms  a  self-enclosed  unit 
contained  in  a  heavy-walled  soft  iron  tube  3.5  inches  in  diameter 
which  serves  both  as  a  return  flux  path  and  a  structural  frame  of  the 
instrument.  To  it  are  bolted  the  two  quadrupole  magnet  systems  that 
guide  the  suspended  dipole  along  the  vertical  axis  and  prevent  its  tip¬ 
ping  over  due  to  the  inherent  instability.  The  complete  seismometer  is 
shown  in  Figure  11  and  some  of  its  essential  parts  are  shown  disassembled 
in  Figure  12.  The  seismometer  fits  in  a  standard  4.75  in.  O.D.  steel 

IV- 7 


3rthur  2D. little. Jfitr. 


rv-8 


IV -9 


case  (C)  which  serves  again  as  a  magnetic  shield.  The  construction 
is  simplified  somewhat  by  the  fact  that  the  instrument  need  not  be 
levelled  to  the  narrow  limit  required  by  the  horizontal  seismometer. 

The  present  model  tolerates  deviations  from  the  vertical  of  approximately 
±1  degree  and  we  are  trying  to  extend  this  to  ±2  degrees  by  strengthen¬ 
ing  the  quadrupole  systems  Q.  Consequently,  we  do  not  plan  to  use  any 
remote  levelling  system  for  the  vertical  seismometer. 

The  electrooptical  readout  system  (EO)  is  mounted  on  the  top  of  the 
instrument.  It  consists  of  a  light  source  and  a  dual  photocell  of  the 
same  type  as  in  the  horizontal  instrument.  The  vertical  motion  of  a 
disc-shaped  light  obstruction  attached  to  the  upper  end  of  the  seismic 
mass  controls  the  distribution  of  light  between  the  two  halves  of  the 
CdS  photodetector. 

Referring  to  Figure  12  we  see  the  suspension  magnet  system  (partially 
opened)  showing  the  upper  polepiece  (S)  attached  to  a  ring-shaped 
Alnico  V  permanent  magnet.  The  circular  polepiece  (N)  represents  only  a 
fraction  of  the  theoretical  hyperboloid  contour;  however,  this  is  found 
to  be  satisfactory  in  practice.  The  upper  and  lower  polepieces  (S)  are 
provided  with  0.25  in.  diameter  axial  holes  to  permit  the  passage  of 
the  seismic  mass  assembly  (M)  through  them.  The  seismic  mass  consists 
of  several  parts  held  together  by  a  length  of  a  0.187  in.  O.D., 

0.005  in.  wall  stainless  steel  tubing  (t).  In  the  center  of  the  tube 
is  fixed  the  suspended  dipole  which  has  a  form  of  cylinder  0.177  in. 

O.D. ,  0.375  in.  long  made  of  sintered  Alnico  VIII;  this  material  has 
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sufficiently  high  coercive  force  (~1600  Oe)  to  minimize  the 
undesirable  effect  of  position-dependent  induced  magnetization. 

The  upper  and  lower  ends  of  the  stainless  steel  tube  hold  two 
graphite  rods  (g),  0.187  in.  O.D.,  which  serve  as  diamagnetic  guides 
for  the  suspended  seismic  mass.  The  mass  is  augmented  by  two  brass 
slugs  (m)  which  fit  closely  in  the  tube  and  the  weight  of  which  is 
adjusted  roughly  so  that  Equation  (12)  is  fulfilled.  Fine  adjustment  of 
the  gradient  is  done  by  means  of  magnetic  shunt  screws  adjacent  to  the 
Alnico  V  ring  magnets  (not  shown  in  the  picture).  Fine  adjustment  of  the 
second-order  term  (Eq.  19)  which  determines  the  spring  constant  of  the 
suspension  is  made  by  means  of  a  thin  coil  located  in  the  plane  of  the 
ring  (N);  by  adjusting  the  current  through  this  coil  the  period  can  be 
set  to  any  desired  value.  At  present  the  period  is  set  to  approximately 
2.5  seconds.  Another  coil  ("motor"  coil)  serves  to  produce  a  small 
variable  gradient  for  the  purpose  of  applying  a  control  force  on  the 
suspended  dipole.  The  complete  set  of  coils  (not  shown  in  Fig.  11)  fits 
coaxially  with  the  seismic  mass  assembly  in  the  working  space  of  the 
suspension  magnet  assembly. 

The  feedback  system  for  the  vertical  seismometer  is  similar  to  that 
for  the  horizontal  instrument  in  principle  but  differs  in  several  details. 
The  circuit  diagram  is  shown  in  Fig.  13.  To  allow  for  a  possible  4 : 1 
change  of  gain  with  period,  a  gain  step  of  16  (24  db)  was  used.  The  gain 
and  the  frequency  of  the  gain  step  are  adjustable.  The  output  current  is 
fed  back  to  the  motor  coil  of  the  seismometer.  The  feedback  loop  includes 
an  integrating  amplifier  which  provides  high  gain  at  very  low  frequencies, 
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below  the  band  of  interest.  The  effect  of  this  integration  is  to  lower 
the  overall  response  at  these  low  frequencies  and  thus  minimize  the 
drift  of  the  center  point  of  suspension. 
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V. 


BOREHOLE  INSTRUMENT  PACKAGE 


For  borehole  installation  the  two  horizontal  seismometers  and  the 
vertical  seismometer  are  combined  in  one  unit  schematically  outlined  in 
Figure  14.  This  package  is  still  under  construction  and  has  not  been 
tested  as  yet.  It  consists  of  a  stainless  steel  tube  5.50  in.  O.D., 

.120  in.  wall,  approximately  40  in.  long  with  a  welded  bottom  and  an 
O-ring  sealed  header  with  watertight  cable  feedthroughs.  The  three 
instruments  are  stacked  in  this  tube  and  bolted  together  with  .2  in. 
spacer  rings.  The  upper  ring  has  an  expandable  section  to  assure  firm 
coupling  to  the  wall  of  the  tube  which  in  turn  is  eventually  coupled  to 
the  casing  of  the  borehole.  We  are  presently  planning  to  use  Ser.  SW 
borehole  casing  of  6  in.  I.D.  and  6-5/8  in.  O.D.  for  down-hole  testing 
of  the  instrument.  The  upper  part  of  the  instrument  package  is  coupled 
to  the  borehole  casing  by  a  suitable  hole  lock.  We  have  not  yet  made 
the  final  decision  whether  to  use  a  simple  tension  release-operated  lock 
or  a  motor-driven  one.  Three  multi-conductor  cables  will  connect  the  in¬ 
strument  with  e  junction  box  and  the  electronics  console  on  the  surface. 
Even  though  we  expect  to  have  a  "dry"-cased  borehole  for  the  tests,  the 
instrument  package  and  the  cabling  will  be  constructed  to  be  watertight. 
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VI.  PRELIMINARY  TEST  RESULTS 


A.  FREQUENCY  RESPONSE  OF  THE  SEISMOMETERS 

Tests  reported  below  are  of  a  preliminary  nature  inasmuch  as  they 
have  been  performed  on  f irst-of-a-kind  experimental  instruments  and 
their  purpose  was  mainly  to  ascertain  whether  or  not  the  instruments 
perform  according  to  the  design  goals.  Also,  the  parameters  of  filter 
networks  that  will  be  eventually  used  for  shaping  of  the  passband  have 
not  been  finalized  at  the  time  of  the  tests.  For  this  reason  the  fre¬ 
quency  response  of  the  horizontal  seismometer  was  tested  without  any 
filters  outside  the  feedback  electronics. 

The  horizontal  instruments  were  tested  on  a  rocking  platform  pro¬ 
ducing  an  acceleration  cf  a  constant  amplitude  (approximately  6.4  x  10  ^  g) 
and  variable  frequency  (from  .008  to  1  Hz).  Alternately,  we  applied  a 
low  frequency  test  voltage  on  the  PZT  crystals  actuating  the  tilt  plat¬ 
form.  High  amplitude  of  the  test  signal  was  necessary  in  order  to  main¬ 
tain  a  sufficiently  high  signal-to-noise  ratio  at  the  high  ground  noise 
level  prevailing  at  the  Acorn  Park  laboratories. 

A  typical  test  result  is  shown  in  Figure  15.  Both  amplitude  and 
phase  response  are  plotted  on  the  diagram.  The  instrument  was  a  hori¬ 
zontal  seismometer  No.  C-2  operated  with  the  feedback  circuit  described 
earlier  (Section  III). 

Results  of  a  similar  test  made  with  a  vertical  seismometer 
No.  VS-2  are  shown  in  Figure  16.  Since  we  did  not  have  a  shake  table 
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PHASE,  DEGREES  RELATIVE  RESPONSE 


FIGURE  15  FREQUENCY  RESPONSE  OF  THE  HORIZONTAL  SEISMOMETER  C-2 

WITH  FEEDBACK,  TESTED  AT  CONSTANT  ACCELERATION 
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FIGURE  16  FREQUENCY  RESPONSE  OF  THE  VERTICAL  SEISMOMETER  VS  2  WITH 
FEEDBACK,  TESTED  AT  CONSTANT  ACCELERATION 


operable  over  the  required  low  frequency  range,  we  used  an  indirect  method. 
A  test  signal  derived  from  a  low  frequency  function  generator  (Hewlett- 
Packard  Mod.  202A)  was  fed  into  one  of  the  motor  coils  to  generate  a 
frequency-variable  force  of  essentially  constant  amplitude  acting  on  the 
seismic  mass.  Consequently,  tests  performed  in  this  way  yield  frequency 
response  versus  acceleration  characteristics.  The  curves  plotted  in 
Figure  16  were  obtained  with  the  vertical  seismometer  operated  with  a 
feedback  circuit  according  to  Figure  13.  A  simple  low  pass  filter  with 
a  cutoff  frequency  of  approximately  1  Hz  was  used  in  this  test  at  the 
signal  output. 

B.  SAMPLE  SEISMOGRAMS 

During  the  laboratory  tests  at  Acorn  Park  (Cambridge,  Massachusetts) 
teleseismic  events  were  occasionally  recorded.  Because  of  the  seis- 
mologically  unfavorable  characteristics  of  the  site  (swampy  ground 
without  connection  to  bedrock,  combined  with  a  very  high  level  of 
natural  and  man-made  ground  noise)  the  instruments  could  never  be  operated 
at  high  gain.  Consequently,  only  large  events  (M>  6.5)  were  recorded. 

Also,  since  no  effort  was  made  to  operate  the  instruments  continually, 
many  events  undoubtedly  escaped  ovr  attention.  Under  these  conditions 
obtaining  a  good  record  was  a  matter  of  chance.  A  sample  of  one  of  the 
best  recordings  obtained  with  the  horizontal  feedback  seismograph  is 
shown  in  Figure  17.  The  event  was  an  earthquake  (M  ft  7)  off  Hokkaido, 
Japan,  which  occurred  at  10:39:01.6  on  16  May  1568.  The  epicenter  was 
identified  from  the  C.G.S,  -  P.D.E.  cards.  The  record  is  good  enough 
to  permit  clear  identification  of  most  of  the  typical  phases  present. 
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FIGURE  17  PART  OF  THE  RECORDING  OF  A  TELESEISMIC  EVENT 

MADE  WITH  THE  C-2  HORIZONTAL  FEEDBACK  SEISMOGRAPH 


APPENDIX  A 


ANALYS IS  OF  FEEDBACK  VERTICAL  SEISMOMETER 


In  the  development  of  the  vertical  seismometer 
using  magnetic  suspension  and  feedback, calcula¬ 
tions  on  the  response  of  the  system  to  the  fol¬ 
lowing  input  signals  have  been  made:  vertical 
displacement  of  the  '.rarae;  tilting  of  the  frame; 
Brownian  motion  of  the  seismic  mass,  and  elec¬ 
tronic  noise. 
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I. 


INTRODUCTION  AND  SUMMARY 


The  attached  series  of  calculations  is  an  attempt  to  collect 
in  one  document  the  analysis  that  has  been  done  over  the  last  year  on 
the  development  of  a  fedback  magnetically  supported  vertical  seismom¬ 
eter.  It  should  be  emphasized  that  this  is  a  preliminary  and  incom¬ 
plete  attempt.  The  present  work  still  requires  editing  to  be  concise 
and  internally  consistent.  In  addition,  the  work  requires  extension  so 
that  hard  and  fast  design  decisions  can  be  made. 

Two  tentative  designs  are  sketched;  one  employing  positive 
feedback  to  increase  the  natural  period  of  the  instrument,  and  negative 
feedback  to  reduce  drift.  This  design,  described  as  design  A  in 
section  XI,  requires  wide  dynamic  range  in  mass  displacement,  and  its 
output  is  sensitive  to  detector  sensitivity  and  amplifier  gain  varia¬ 
tions. 

The  second  design,  suggested  by  Martin  Cohen,  is  fed  back 
negatively  over  the  full  range  of  operating  frequencies.  Its  gain  is  a 
function  only  of  the  force  constant  of  the  instrument  and  is  independent 
of  detector  sensitivity  and  amplifier  gain. 
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II.  FEDBACK  VERTICAL  SEISMOMETER 


FIGURE  II-l.  DIAGRAM  OF  FEDBACK  VERTICAL  SEISMOMETER 

where 

Z^(t)  ■  the  time  varying  displacement  of  the  frame  (cm) 

Z(t)  *  the  time  varying  displacement  of  the  mass  m 

with  respect  to  the  frame  (cm) 

Zo  -  the  equilibrium  position  of  the  mass  m  with 

respect  to  the  frame  (cm) 

X(t)  -  Z(t)  -  7o  (cm) 

j3  -  detector  sensitivity 

(volts/cm  mass  displacement) 
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n  ( t ) 


A 

0(t) 


B 


F 

o 


mg 


■»  equivalent  input  electrical  noise  to  the 
system  volts 

■  amplifier  gain 

■  output  signal  volts 

=  feedback  factor  cm  ^/volt 
a  feedback  force  (dynes) 

■  static  supporting  force  (dynes) 

=  static  gravitational  force  (dynes) 


-mg02 (t j 
2 


0(0 

F(0) 

e 


tilt  force  (dynes) 

tilt  angle  relative  to  the  vertical  (radians) 
Brownian  force  (dynes) 
temperature  (degrees  absolute) 


Considering  the  total  vertical  acceleration  of  the  seismic  mass  m  to  be 
the  acceleration  of  the  frame  plus  the  acceleration  of  the  mass  relative 
to  the  frame,  one  can  write  a  general  force  balance  equation  of  the  form 

mX  +  Q*  -  -mZ,  +  F  +  F(0)  +  F_(t)  -mg  (1  -  ^  )  (II-l) 

to  D  L 
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III.  MAGNETIC  F II7, LD  DESIGN  AND  ADJUSTMENT 


If  the  magnetic  field  is  configured  so  that 

H(X)  =  Hq  £{a  +  b(t)}  X  -  |  X2J  (III-l) 

and  further,  if  the  seismic  mass  m  is  given  a  magnetic  dipole  strength 
of  M>,  then  the  vertical  magnetic  force  on  the  mass  will  be  given  by 


F  +  F  (t)  -  -  M  [a  +  b(t)  -  c  .  X(t) ]  (III-2) 

o  B  dX  o 


where 

F  ■  static  supporting  force  (dynes) 
o 

F  (t)  «*  time  varying  fedback  force  (dynes) 

B 

U  ■  magnetic  dipole  strength  of  seismic  mass 
(dyne -cm  .  gauss"l) 

H  «*  magnetic  field  at  equilibrium  point, 
i .e .  ,  X**o  (gauss) 

a  •»  static  field  strength  adjustable  parameter 
(cra"l) 

b(t)  ■  dynamic  field  strength  adjustable  parameter 
(cm"l) 

c  «  static  field  curvature  adjustable  parameter 

(cm" 2) 


Referring  to  equation  (II-l)  and  equating  static  terms,  at  X=o,  one 


obtains 

F  -  H  H  a  =  mg 
o  o 


(III-3) 
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To  adjust  the  operating  point  at  X*o,  one  adjusts  the  static  parameter  a 


so  that 


PR 

|IH 


(III-4) 


Assuming  this  adjustment,  equation  (II-l)  can  be  rewritten  in  the  form 


m  m 


.  i  +  m  +  4a  +  3n(t) 


(III-5) 
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IV.  RESPONSE  TO  SINUSOIDAL  FRAME  DISPLACEMENT 


In  the  absence  of  extraneous  signals 


b(t)  -  X(t)  •  Bab 


(IV-1) 


Therefore,  from  equation  (III-5) 


M>H 

X  +  —  X  +  — — ■  (c  -  &AB)X  =  Z, 
mm  r 


(IV-2) 


Taking  the  transform  of  equation  (IV-2)  and  solving  for  the  transfer  gain 
of  the  system,  one  obtains 


defining 


0(a)  _ 
Z(s) 


^  +  2  .  *J!°L  (1 . 

mm  c 


j¥|; 


(undamped  resonant  frequency 
of  the  unfedback  system) 


(damping  factor  for  unfedback 
system) 


(1  -  feedback  parameter 


0(8) 

Zf(s) 


2  2  2 
s  +26ujs+y  <ju 
o  o 


(IV-3) 


(IV -4) 
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(Db  =  20  log 


letting  s  *»  icu  one  obtains  an  expression  for  the  frequency  responie  of 


the  system 


(IV-5) 


This  function  is  plotted  in  Figure  IV-1  for  an  assumed  value  of  6  =  0.1 
and  0.1  <  y  <2.0 


FIGURE  IV-1  FEDBACK  SEISMOMETER  RESPONSE  TO  SINUSOIDAL 
DISPLACEMENT  SIGNALS 
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Referring  to  the  definition  of  y,  the  feedback  factor,  it  is  seen  that 
when  y  ■  1,  there  is  no  feedback  and  the  system  has  the  characteristics 
of  a  simple  pendulum  with  damping  factor  6  and  undamped  resonance  fre¬ 
quency  •  Phasing  the  amplifier  so  that  it  reduces  y  reduces  the 
undamped  resonance  frequency  of  the  system  and  increases  the  apparent 
damping  factor.  There  is  a  limit  on  loop  gain  in  that  when  y  =*  o,  the 
system  becomes  exponential  rather  than  oscillatory,  and  if  the  gain  is 
further  increased,  the  system  becomes  exponentially  diverge 

Phasing  the  feedback  so  that  y  is  increased  beyond  unity 
corresponds  to  negative  feedback  in  that  the  system  becomes  stiffer, 
damping  is  reduced,  and  undamped  resonant  frequency  is  increase  d. 
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V.  BROWNIAN  MOTION  ANALYSIS 


Referring  to  equation  (III-5)  and  under  the  assumption  that 
b(t)  =  X(t)  0AB,  then  with  a  Brownian  force  F(0),  one  can  write 


M.H 

i i  +  -  x  +  — —  (c  -  pab)  x  =  F^e,t^ 

mm  m 


(V-l) 


F  (0)  , t)  /m _ 

52  +  a  s  +  ^  (1  .  M  } 

m  m  c 


(V-2) 


f (0.s)/m _ 

2  2  2 
3  +26ojs+Yu) 
o  o 


(V-3) 


PAX(s) 


2  2  2 
s  +26ujS+Yoj 
o  o 


(V-4) 


0(0, cu) 


-2i6uj  w  +  Y^ou^ 
o  o 


(V-5) 


0(0,  uj)  = 


pAf  (0 ,  uj) 
2  2 

Y  id  m 
o 


) "¥•(*) 


(V-6a) 


It  will  be  shown  that 


f(9,w)  w  {4cyk0}‘ 


where  Of  is  the  damping  coefficient 

k  is  the  Boltzmann  constant 

0  is  the  absolute  temperature 


and  that 


3rthur  ZH.littlr.Jnr. 


) 


V 


0(8uj) 


^[40-ke}1 
2  2 

V  id  m 
o 


(V-6b) 


This  output  is  plotted  as  a  function  of  frequency,  id,  in  Figure  V-l. 

From  Figure  V-l  it  is  obvious  that  as  positive  feedback  is  increased, 
Brownian  noise  voltage  output  is  increased.  It  is  not,  however,  obvious 
that  one  is  loosing  in  signal-to-noise  ratio.  Signal-to-Brownian  noise 
ratio  can  be  calculated  by  dividing  equation  (IV-5)  by  equation  (V-6) 


one  obtains 


-U) 


2  2  Zf  <“> 


0A 


V  UL, 


s  _  0(CD}. 

N  0(0, (d) 


(■  •*>'?(*) 


BAf(0,ip) 

2  2 
my  <d 


-mu/ Zf  (id) 

[2ake]^ 


(V-7) 


From  this  relationship,  one  concludes  that  the  signal -to-Brownian  noise 
ratio  is  unaffected  by  the  adjustment  of  the  feedback  parameter.  It  is 
a  function  of  frequency  (id),  temperature  (0),  the  unfedback  mechanical 
damping  constant  a,  and  input  signal  spectrum  level  Z(u>). 

The  seismometer  designer  can  control  and  maximize  signal-to- 
Brownian  noise  ratio  by  minimizing  a  and  0  and  maximizing  m. 
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RELATIVE  noise  VOLTAGE  output  (Db) 


VI.  DERIVATION  OF  THE  BROWNIAN  FORCE  SPECTRAL  DENSITY 


Consider  a  simple  spring  mass  system  characterized  by  the 


constants  m,  a  and  k  so  that 


my  +  ay  +  ky  ■  B(t) 


(VI-1) 


Taking  the  transform  of  this  equation,  one  obtains  an  expression 
for  the  transform  of  the  displacement  of  the  mass  m  in  terms  of  the  trans¬ 
form  of  the  applied  force,  that  is: 


b(s) 


(VI-2) 


2  .  as  k 
s  +  —  +  ~ 
m  m 


Letting  s  =*  iu>  and  solving  for  the  squared  magnitude  of  the  quantity  on 
the  right,  one  obtains  an  expression  for  the  power  spectral  density  cf 
|y2(f)j  under  the  condition  that  B(T)  is  a  random,  stationary  ergodic 
function^  that  is 


d  2 
df  y 


(VI-3) 


3  to12'*’'1* 

2TTm  < 


/k  2  ,  io\/k  2  icA 

\m  m  f  ym  m  j 


(1)  Measurement  and  Analysis  of  Random  Data,  Bendat  and  Piersol,  John 
Wiley  &  Sons,  New  York,  N.  Y.  1966.  Chapter  1. 
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Formally  integrating  equation  (3),  one  obtains 


Since  .fL  b  (uj)  must  be  a  positive  definite  real  quantity 
du) 


(VI-6) 
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2 

Since  b  (oj)  is  a  positive  real  quantity 


However , 


Therefore, 


^  k  /  -  h  k  0 

where  k  is  the  Boltzmann  constant 

6  is  the  absolute  temperature 


(VI-10) 


(VI-11) 


(VI-12) 
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This  implies  that  d/dou  is  a  function  of  ot  and  0  only  independent  of 
k,  m  or  u).  This  then  states  that  the  Brownian  force  power  density 
spectrum  is  a  constant  for  fixed  damping  coefficient,  a  and  for  fixed 
temperature  0  and  this  is  given  by 


(f)  »  AakB 


(VI-13) 


and 


b(f) 


Uok©} 


h 


(VI-14) 
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VII.  TILT  SIGNALS 


Referring  to  equation  (III-5)  a  generalized  force  balance 
equation  can  be  written  in  the  form 


i'  +  sa  +  i-a-x  .  +  S>  b(t) 

mm  m  m 

(VII-1) 

for  tilt  signals, 

F(t)  -  E&  0 (t) 2 

(VII-2) 

given  that 

b(t)  -  X(t)  0AB 

(VII-3) 

0(t)  -  0A  X(t) 

(VII-4) 

(62  +  26u>  6  +  y2  iu2)X(6)  =  f(6) 

(VII-5) 

(VII-6) 

(6  +  26  ui  6  +  y  u>  ) 

o  o' 


0(w) 


gA  g0  (uu/2) 

„  2  2 
2y  u> 


(VII-7) 


This  equation  is  of  the  same  form  as  equation  (VI-6a)  and  will  therefore 
lead  to  the  same  conclusions,  i.e.,  signal-to-tilt  noise  ratios  are 
unaffected  by  feedback. 

Tilt  Frequency  Doubling 

An  interesting  aspect  of  the  instrument's  tilt  sensitivity  is 
the  frequency  doubling  plus  dc  offset  that  occurs  with  a  periodic  tilt 
signal . 
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If  one  assumes  a  tilt  of  the  form 


0(t)  -  tQ  sin  —  (VII-8) 

then 


„  .  .  mg  .2.  .  mg  .2  .  2  u>t 

F (t)  =  2*  t  (t)  -  -f-  0O  sin  T 

(VII-9) 

F (t)  =  ^  t>?Q  {l  -  cos  cut } 

(VII-10) 

Therefore,  for 
that  has  a  dc 
term  is  given 


a  periodic  tilt  at  frequency  u>/2  one  gets  a  forcing  function 
component  plus  a  periodic  component  at  frequency  u>.  The  dc 

by  2 

PA  g  t 

0(o)  =  - 5—7  (VII-11) 


2  2 
2y  uj 

o 
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VIII.  ELECTRONIC  NOISE 


Referring  to  Figure  II-l,  it  can  be  seen  that 
0(t)  -  (BX(t)  +  T](t)  )A 
0(s)  -  {PX(s)  +  T](b)  }A 


b(t)  -  BO(l) 


F,  -  b(t)  =  (iH  B  0(t) 
bo  o 


„  HH  cX  HH  BA 

X  +  ^  +  — T—  “  (0X(t)  +  T](t)) 


in  ra 


m 


+  ^  (>-¥)} 


M-H  BA 
o 

m 


T)<S> 


HH  BA 

D(s)  -° 


m 


X(s)  - 


2  O' 

S  +  —  s  + 

ra 


Qjsi 

Tl(s) 


RP-**) 


g2+-g+ 

ra  m 
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(VIII-1) 

(VIII-2) 

(VIII-3) 

(VIII-4) 

(VIII-5) 

(VIII-6) 

(VIII-7) 


3rthur  2D.littlf.llnr. 


oM 

T'((0) 


(VIII-8) 


( x 88  *  l) A 
2  2  2 
-u)  +  216  a)  id  +  V  a) 
o  o 

Here  again,  one  sees  the  output  voltage  due  to  noise  input  increases 
as  Y  “*  0  (i.e.,  positive  feedback):  however,  the  signal-to-noise  ratio 
will  again  be  independent  of  che  amount  and  type  of  feedback  used. 
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IX.  ROOT  LOCUS  AND  TRANSIENT  RESPONSE 


Referring  to  equation  (IV-3) 


O^s) 

Zf(s) 


PA 


s 


2  a  2  k 
s  +  -  s  +  v  - 
in  m 


(IX- 1] 


a  2  k 

Upon  substitution  of  6  ■  „ .  ,  u)  •»  —  ,  one  obtains 

v  o  2m  uu  o  m 

o 


0^8) 

z^T 


PA 


s 


2  2  2 
s  +26  ujs+v  u) 
o  o  o 


(IX-2 


The  denominator  can  be  factored  in  that 


±{ 


or  ,  /  .  *2  2  .  2  2  1 

-26  uj  +146  u)  -4yu)  ( 


o  o  —  l  o  o 


s  =  cu 


[■6„± 


01(s) 

Zf(s) 


*  PA 


_s+Vo  +  {6o  ’  y2}^]  8  +  6o  "  {6o  ■  y2}  ] 


The  function  in  (IX-3)  will  have  poles  at 


S  ■  UJ 


There  is  a  series  of  points  of  interest: 
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real  and  >  6. 


CASE  I 


2 

Y  <  0 


.2  2  £ 
6  -  y 

o 


The  solution  is  therefore  made  up  of  exponen¬ 
tials  of  the  form 


■[26  +  c6  ]t 

o 


+e6  t 

0,0 

and  e 


CASE  II 


The  system  therefore  is  unstable,  due  to  the 
positive  exponential  term. 

2 

Y  =0  Under  this  condition  the  function 
has  a  pole  at  the  origin  and  at  -26q,  solu¬ 
tions  will  contain  a  single  solution  of  the 

form  -26  t 

o 

e 

The  system  is  stable  and  approaches  equilib¬ 
rium  exponentially. 


CASE  III 


2  2 
0  <  Y  <  6 

o 


1 .2  2 

'6  -  y 

o 


real  and  <  6 


The  solution  is  therefore  made  up  of  exponen¬ 
tials  of  the  form 

-[26  -  e6  ]t  -e  6  t 
0  0  o 

e  and  e 

The  system  is  therefore  stable  in  that  it 
approaches  equilibrium  exponentially. 
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CASE  IV 


2  2 

Yq  ■  6q  ,  the  system  has  a  pole  of  order  2  at 

-6  and  solutions  will  be  of  the  form 
o 

-6  t 

1  -  t  e  ° 

The  system  is  stable  but  now  approaches  equilib¬ 
rium  after  one  zero  crossing. 


CASE  V 


Y  >  6 
o  o 


The  system  has  poles  at 

ft 


s  ■  -<u  6  +  itu  <  Y^  ~  6  i 

o  o  -  O  l  O  / 


Solutions  will  be  of  the  form 


-Y  t 


e  °  sin  (Uo  |y2  -  6^  |  t 


The  system  is  stable  and  will  approach  equilib¬ 
rium  with  a  damped  sine  wave  response 

The  root  locus  can  now  be  sketched  in  the  complex  s  plane  as  shown: 


FIGURE  IX- 1  ROOT'  LOCUS  FOR  TRANSFER  GAIN 
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FIGURE  IX-2  FEDBACK  SEISMOMETER  TRANSIENT  RESPONSE 
FOR  STEP  MOVEMENT  OF  FRAME 
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X. _ r  IIFT 

Substituting  equation  (III-2)  into  equation  (II-l)  in  the 
absence  of  ground  motion,  Brownian,  tilt  and  electronic  noise,  one  obtains 


mx  +  ox  “  M-H  Ca  +  b(t)  -  c  •  x(t)3  -  mg 
o 


Now  let  us  assume  that  Hq,  a  or  c  are  slowly  varying  functions  of  time, 
then  as  in  •*  0 


a  +  b(t)  -  c  X(t) 


(X-l) 


However,  from  equation  IV-1  b(t)  =  PAB  X(t) 


a  +  (PAB  -  c)  X(t) 

x(t) 


i2S—  -  a  i  — - 

M>Hq  j  3AB  -  c 


(X-2) 

(X-3) 


To  assure  that  X(t)  ■*  0  at  dc  (ur  -*  o)  it  is  necessary  to  have  PAB  "*  oo 
as  w  "  o,  This  implies  very  strong  negative  feedback  at  low  frequencies. 

If  we  compare  the  instrument  drift  with  and  without  feedback, 


we  have 


X(t)  feedback 
X(t)  open  loop 


uu 


o 


P 

|^Ho 


1 

c  -  PAB 

“T 

c 


(X.4) 
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X(t)  feedback 

X(t)  open  loop  ^  _  Q 


1  =  1 
1  .  M®.  °  \  O) 


(X”5) 


The  drift  of  the  instrument  is  therefore  reduced  by  a  factor  of 


V  -*  o) 
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XI.  FEEDBACK  CIRCUIT  DESIGN  A 


Consider  a  circuit  of  the  type  shown  in  Figure  XI-1. 


FIGURE  XI-1.  STABILIZED  FEDBACK  VERTICAL  SEISMOMETER 


In  the  circuit  sketched  one  feeds  back  negatively  at  low  fre¬ 
quencies  and  positively  at  high  frequencies.  A  sketch  of  the  Bode 
diagram  for  this  system  is  shown  in  Figure  XI-2. 
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FIGURE  XI-2.  BODE  DIAGRAM  FOR  STABILIZED  SEISMOMETER 
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The  upper  diagram  in  Figure  XI-2  corresponds  to  the  case 

where  is  a  simple  low-pass  filter  with  a  corner  frequency  at 

2  x  10  u)  ,  and  B.  .  ,  is  a  high  pass  filter  with  a  corner  frequency 
o  nigh 

-2  . 

at  5  x  10  idq.  The  lower  diagram  indicates  the  change  in  response 
if  the  filters  cut  off  at  the  rate  of  40  db/decade  rather  than  20  db/ 
decade.  The  diagrams  are  for  the  overall  input-output  frequency 
response. 

It  should  be  emphasized  that  in  this  design  the  overall  fed- 
back  dynamic  response  of  the  mass  is  equivalent  to  that  of  a  critically 
damped  simple  pendulum  with  undamped  resonance  frequency  equal  to 
0.2  times  the  undamped  resonance  frequency  of  the  unfedback  system. 

In  this  design,  the  mass  moves  with  respect  to  the  frame. 
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XII.  FEEDBACK  CIRCUIT  DESIGN  B 


Because  of  dynamic  range,  linearity,  temperature,  stability 
and  gain  constancy  considerations,  Martin  Cohen  suggests  that  the 
system  should  be  negatively  fedback  over  the  entire  operating  range  so 
that  the  seismic  mass  is  constrained  not  to  move  over  the  frequency 
range  of  interest  and  further,  that  only  a  portion  of  the  shaping  be 
done  within  the  loop  with  final  shaping  and  gain  adjustment  being 
accomplished  by  post-amplification  and  filtering.  A  design  of  this 
type  is  sketched  in  Figure  XII- 1. 


FIGURE  XII-1.  NEGATIVELY  FEDBACK  VERTICAL  SEISMOMETER 
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FIGURE  XII-2.  ACCELERATION  RESPONSE  OF  GAIN  COMPENSATED  NEGATIVELY 

FEDBACK  SEISMOMETER 


The  device  sketched  will  have  an  acceleration  response  of  the  form 


oiiL 


s*  zf(s) 


BA(s) 


2  O' 

s  +  -  s  + 

m  m 


o_l  ^  .  PA(s)B(s).j 


(XII-1) 


If  in  the  frequency  range  of  interest  A(s)  ■  A  {s  +  a^,  then 


BA  {s  +  a} 
o 


0(8)  .  _ 

2  .  .  0  /  BA  B(s)  c  \  M-H  c  /  BA  B  \ 

8  zf(s>  2 - a- »)+ -a- (l  + -a- . 

\m  a  m  /mV  c  ) 


(XII-2) 


0(s) 


BA  {s  +  a} 
o 


s  Zf(s) 


2 

s  + 


/  Y^O)  ?  -  1  \  2  2 

( 2  6uu  +  - 2 - )  s  +  yu 

\  o  a  /  o 


(XII-3) 


Since  the  constant  a  can  be  chosen,  the  damping  at  the  fedback 
resonance  can  be  adjusted  to  obtain  a  critically  damped  system.  Tht  Bode 
diagram  for  this  system  is  shown  in  Figure  XII-2.  As  seen  from  this  dia¬ 
gram,  the  acceleration  response  over  the  range  of  frequencies  of  interest 
is  flat. 

The  transfer  gain  for  acceleration  in  this  region  is  given  by 

0(s)  BAa 

2  2  2 

S  6f  (s)  YV 

and  therefore  a  function  of  the  force  constant  of  the  feedback  only. 


BAa 


HH  c  /  BA  Ba\ 


HH 


(XII-4) 
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